Model-independent WIMP Characterisation using ISR 
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The prospects of measuring the parameters of WIMP dark matter in a model indepen- 
dent way at the International Linear Collider are investigated. The signal under study 
is direct WIMP pair production with associated initial state radiation e + e~ — > XX7- 
The analysis is performed in full simulation of the ILD detector concept. With an 
integrated luminosity of C = 500 fb _1 and realistic beam polarizations the helicity 
structure of the WIMP couplings to electrons can be determined, and the masses and 
cross sections can be measured to the percent level. The systematic uncertainties are 
dominated by the polarization measurement and the luminosity spectrum. 

<j_) _ 1 Radiative WIMP production in e + e collisions 



New Weakly Interacting Massive Particles (WIMPs) with masses in the order of M x ~ 
100 GeV are predicted by several extensions to the SM of particle physics. With their weak 
strength interactions, these particles are natural candidates for the observed abundance of 
cosmological Dark Matter (DM). 

If these particles were produced at colliders, they would leave the detector invisibly 
without any further interaction. Their parameters could be inferred indirectly by the analysis 
of cascade decays if other new particles exist in the kincmatically accessible mass range. 
Alternatively the direct pair production of WIMPs with associated initial state radiation 
e + e~ — > XXI can be employed to determine the WIMP properties from the observed photon 
spectrum. It has to be noted, that even if the detection via cascade decays is possible, 
the single photon plus missing energy signature provides an additional measurement of the 
WIMP candidate. 

The rate of radiative WIMP production at an e + e~ collider can be estimated model- 
independent without any assumptions on the dynamics of the interaction involved [T] . With 
only one new stable particle responsible for the observed DM content in the universe, the 
production cross section for WIMP pairs with associated ISR can be written in the limit of 
non-relativistic final state WIMPs as: 

da _ aK e a an 1 + (1 - xf 2J a / 4Mg \ 1 +h 

dxdcosO~ 16tt xsin9 2 [ x + ' [ (I - x)s J ' [ ) 

with the candidate mass M x , the candidate spin S x and the center-of-mass energy squared s. 
The double differential cross section is expressed in the dimensionless variables x = and 
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of the emitted photon. In Equation [1] Jo is the quantum number of the dominant partial 
wave in the production process. In the following, we will refer to the cases of Jq = (s- 
wave) and Jq = 1 (p-wave) production only. The quantitiy n e is the "annihilation fraction" 
of WIMPs into electrons. The parameter a an provides the overall scale of the production 
cross section and can be inferred from observation when the WIMP is identified with the 
cosmological Dark Matter pQ. 

The annihilation fraction n e implicitcly depends on the helicity of the initial state elec- 
trons. For our analysis wc investigated the following three coupling scenarios [2]: 

• " Equal" : The WIMP couplings are independent of the helicity of the incoming elec- 
trons and positrons, i.e. «(e^,ej) = «(e^, ej) = K(e£,eJ) = /c(e£,e^). 

• "Helicity": The couplings conserve helicity and parity, 

• "Anti-SM": This scenario is a "best case" scenario, since the WIMPs couple only to 
right-handed electrons and left-handed positrons: Ac(e^,e£). 

This "single photon plus missing energy signature" has Standard Model (SM) processes 
with large cross sections for background, the dominant one being radiative neutrino produc- 
tion e + e~ — > uw), which proceeds for high ccnter-of-mass energies primarily via t-channel 
W exchange and hence is strongly polarisation dependent. In addition other SM processes 
like radiative Bhabha scattering e + e~ — > e + e~j or multi-photon final states can mimic the 
WIMP production signature when the accompanying electrons or photons leave the detector 
through the beam pipe, or are not properly reconstructed. Polarised beams can be used to 
significantly reduce these backgrounds and increase the S/B ratio [3J, but the still large 
abundance of background events requires a consideration of sytematic uncertainties from 
the detector measurement and beam parameters. 

2 Data samples and event selection 

To cover a broad range of parameters in terms of candidate masses, partial waves and 
coupling structures, only the SM background has been generated and simulated explicitly. 
The signal contribution to the data is obtained by reweighting the irreducible SM 1/1/7 
background, which is indistinguishable from the signal on an event- by- event basis, with 
the cross section ratio w S i g (E 7 ) = do ^xx7) j da ^' yS> of radiative WIMP and neutrino pair 
production. The weights are evaluated in terms of the MC photon energy. 

The measurement of the WIMP parameters requires a precise prediction of the SM 
background photon distributions. The expectation is generated by a parametrization of 
an independent background subsample, by succesively correcting the SM prediction for the 
detector energy resolution and selection efficiencies. Remaining differences to the simulated 
detector output from the beam energy spectrum and unaccounted detector and reconstruc- 
tion effects are parametrized with a higher order polynominal. From the parametrization of 
the vvy spectra, the signal prediction is generated [2]. 

2.1 Simulation and reconstruction 

The background events have been generated using Whizard j4|. The beam energy spec- 
trum was provided by Guinea Pig [5] for the nominal RDR baseline parameter set for a 
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y/s = 500 GeV machine [BJ. The simulation of the detector response was done for the ILD 
detector concept [7] with the Geant4 [5] based simulation software Mokka [3]. Three 
statistically independent event samples of ~ 50 fb" 1 each have been generated and sim- 
ulated for the dominant irreducible 2/2/7 background. The first two samples are used for 
the background and signal contribution to the data, while the third sample serves as basis 
for the spectrum parametrization. In addition data samples with multi-photon final states 
e + e~ —J- 77(^)7 and the Bhabha background e + e~ — > e + e~7 have been simulated. After 
event reconstruction using the Pandora PFA particle flow algorithm [10] two corrections to 
the reconstructed events have been applied. The reconstruction of high energy photons often 
results in several additional detections of lower energy photons, because the clustering stage 
of the Pandora PFA algorithm tends to split large energy depositions into several smaller 
distrinct electromagnetic clusters. To counter the effect of the cluster fracturing, photon 
candidates are merged with a cone based method. Second, the photon candidate energies 
have been recalibrated with a calibration function accounting for polar angle dependent 
fractional energy losses caused by the segmentation of the ILD calorimeter system. 

2.2 Event selection 

An event is considered signal-like, if it contains at least one high px photon with an energy 
between 10 GeV < E 1 < 220 GeV, and a polar angle constrained to | cos0| < 0.98. The 
condition on the photon energy reduces the abundant low energy ISR from the SM back- 
ground and excludes the massless neutrino final state on the radiative Z return at photon 
energies of E 1 m 241 GeV. For further event selection, additional constraints are set to 
deal with the dominant reducible SM backgrounds, especially radiative Bhabha scattering 
and multi-photon final states. To exclude hadronic and leptonic final states, the maximal 
exclusive energy E y — E V i s , i.e. the full visible energy excluding the selected photon, is con- 
strained to 20 GeV. For further reduction of hadronic final states and Bhabha events, the 
maximal allowed transverse track momentum is pt < 3 GeV. The track momenta can not 
be constrained stronger, as the event selection has to allow for tracks of e + e~ pairs from 
the beamstrahlung background and for track overlays of multi-peripheral 77 — > hadrons 
events (collectively called 77 processes). In Figure QJa) the momentum distribution of tracks 
per bunch crossing from the 77 and beamstrahlung background are shown. On average 0.7 
tracks from 77 processes and 1.5 tracks from beamstrahlung are expected in each event. 
The distributions show a strong peak at low pt determined from the minimal momentum 
for the tracks to reach the tracking region of ILD and fall of rapidly. Only 0.2% of tracks 
have a transverse momentum above 3 GeV. The selection efficieny of the 2/2/7 background 
is on average above 85%, with higher efficiencies in the low energy part of the photon spec- 
trum, where most of the WIMP production signal is located. This higher efficiency for low 
energies translates to a signal selection efficiency well above 90% over the full mass range, 
see Figure QJb). 

2.3 Systematic uncertainties 

While the WIMP mass is predominantly determined from the polarisation independent 
threshold in the photon energy spectrum, the cross section determination requires a pre- 
cise knowledge of the normalisation of the polarisation dependent 1/2/7 background. The 
polarisation measurement precision is assumed to 8P/P = 0.25% [TT] . The precision on 
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Figure 1: (a) Transverse momentum distribution of track overlays from the beamstrahlungs 
background (red/gray) and from 77 processes (black), (b) Mass dependent signal selection 
efficiencies for s-wave (red/gray) and p-wave (black) WIMP production. 



the luminosity is given by the RDR [6] with 1(T 4 . A further source of systematic error is 
given by the knowledge of the beam energy spectrum. Here, its influence is conscrvatcly 
estimated from the impact on the signal spectrum for two different parameter sets, namely 
the nominal RDR and SB-2009 [12] sets. For the cross section measurement which is not 
sensitive to the partial wave quantum number, the selection efficiency obtains an additional 
contribution of uncertainty from the difference between the selection efficiencies of s- and 
p-wave production. 



3 Results 

For the analysis a typical running scenario of the ILC is assumed, where an integrated lu- 
minosity of C = 500 fb" 1 is distributed to four polarisation states with (+|P e - 1; — \P e + 1), 
{-\P e -\;+\P e+ \), (+\P e -\;+\P e+ \) and {-\P e - 1; -\P e+ 1). The odd (equal) sign configura- 
tions obtain 40% (10%) of the delivered luminosity each. The absolute polarization values 
are assumed to be \P e -\ = 0.8 and \P e +\ = 0.3 or \P e +\ = 0.6, respectively. The total 
unpolarised signal cross section in the signal region is set to do = 100 fb throughout. For 
the determination of the cross section and coupling structure the candidate mass is fixed to 
150 GeV. 



3.1 Helicity structure and cross section 

With four different longitudinal polarization configurations, the fully polarized cross sections 
(J {L,R}i ancl hence the helicity structure of the WIMP couplings to the beam electrons, can 
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be determined from the cross section deconstruction [3] 
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0.00 ±0.37 (0.29) 
1.95 ±0.38 (0.25) 








1.99 ±0.10 
0.00 ±0.23 
0.00 ±0.23 
1.95 ±0.29 


(0.08) 
(0.14) 
(0.15) 
(0.16) 



Table 1: Fully polarized cross sections <J{r.l} measured within the "Helicity" WIMP 
scenario and for two different absolute polarizations of electrons and positrons. The quoted 
uncertainties are the squared sum of statistical errors and systematic uncertainties, with the 
bracketed values corresponding to an increased precision on the polarization measurement 
of SP/P = 0.1%. 
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Figure 2: Helicity structure of WIMP couplings in terms of the fully polarised cross sec- 
tions (T{r,l} m ( a ) the "Helicity" scenario and (b) the "Equal" scenario. The values are 
normalised to the unpolarised cross section <jq. 

The results for the "Helicity" coupling scenario are listed in terms of the fully polar- 
ized cross sections ct{_r,l} m Tableland shown in Figure [2ja). Figure HJb) depicts the 
mesurement in the "Equal" scenario. Given the input cross section of 100 fb, the fully 
polarized cross sections can be determined to 20 fb to 40 fb with a positron polarization of 
30%. The uncertainties are reduced to 10 to 30 fb with an increased positron polarization 
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of 60%. This is in particular the case for <jsl,r}, wich is primarily determined by the mea- 
surement (7-| with positive electron and negative positron polarization. In that case the 

SM background is maximally reduced. The dominant source of systematic uncertaity stems 
from the polarization measurement. 

Assuming uncorrclated polarisation errors, a combination of the four individual measure- 
ments provides a determination of the unpolarised cross section ao to a precision of 3 to 5% 
for an positron polarisation of P = 30%, depending on the coupling scenario, see Table [2J 
With an increased positron polarisation of P = 60%, the relative precision is increased to 
2.5% . Again the dominant source of systematic uncertainty comes from the polarization 
measurement. 
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Assumed polarization uncertainty 5P/P = 0.25% 


"Equal" 99.0 ± 2.8 ± 4.3 (± 5.1) 
"Helicity" 99.1 ± 2.3 ± 4.0 (± 4.6) 
"Anti-SM" 99.8 ± 1.4 ± 2.8 (± 3.2) 
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Table 2: Measured unpolarized cross section ao by a combination of cross section measure- 
ments with polarized beams for an integrated luminosity of C = 500 fb -1 . 



3.2 Mass measurement and partial wave 

The candidate mass is measured by fitting template spectra for s+b to the measured data 
spectrum. The mass is determined by the x 2 of the measurement. Depending on the 
coupling scenario, polarization configuration, the candidate mass can be determined to a 
level of < 2%, see Table [3] 

Mass WIMP mass: ± stat. ± SE (sys.) ± S£ (sys.) (total) [GeV] 

[GcV] (P e _; P e+ ) = (0.8; 0.0) (P e _; P e+ ) = (0.8; -0.3) (P e _; P e+ ) = (0.8; -0.6) 

"Helicity" scenario 

120 2.67 ± 0.07 ± 1.91 (3.29) 1.92 ± 0.07 ± 1.89 (2.70) 1.53 ± 0.07 ± 1.89 (2.43) 

150 2.11 ± 0.05 ± 1.47 (2.57) 1.62 ± 0.05 ± 1.46 (2.18) 1.23 ± 0.05 ± 1.45 (1.90) 

180 1.78 ± 0.03 ± 1.00 (2.04) 1.36 ± 0.03 ± 1.00 (1.69) 0.94 ± 0.03 ± 1.00 (1.37) 

210 0.78 ±0.02 ± 0.54 (0.95) 0.67 ± 0.02 ± 0.54 (0.87) 0.59 ± 0.02 ± 0.54 (0.80 ) 

Table 3: Statistical and systematic uncertainties on the measured WIMP masses for an in- 
tegrated luminosity of C = 500 fb -1 in the "Helicity" coupling scenarios for three different 
polarization configurations. 

The dominant source of systematic uncertainty come from the calibration of the overall 
energy scale (5E) and the measurement of the beam energy spectrum (SC). 

An indication of the dominant partial wave in the production process is obtained from 
the x 2 value of fitting s- and p-wave template spectra against the data spectrum. In all 
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studied scenarios the fit converges better for the correct partial wave assumption. However, 
with unpolarised beams, the distinction is not very clear over the full mass range. Utilising 
the possibility of polarised positrons allows to clearly seperate the s- and p-wave production. 
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